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ABSTRACT 

Deep HS T-WFPC2 observations of the stellar population in the glob- 
ular cluster u Cen (NGC 5139) have been used to measure the luminosity 
function of main sequence stars at the low-mass end. Two fields have 
been investigated, located respectively ~ 4.' 6 and ~ T away from the 
cluster center, or near the half-light radius of this cluster (r^ ~ 4.' 8). 
*5 ' The color-magnitude diagrams derived from these data show the clus- 

ter main sequence extending to the detection limit at I ~ 24. Infor- 
+3 '. mation on both color and magnitude is used to build the luminosity 

functions of main sequence stars in these fields and the two independent 
determinations are found to agree very well with each other within the 
observational uncertainty. Both functions show a peak in the stellar dis- 
tribution around Mj ~ 9 followed by a drop at fainter magnitudes well 
before photometric incompleteness becomes significant, as is typical of 
other globular clusters observed with the HST. This result is at variance 
with previous claims that the luminosity function of u Cen stays flat at 
low masses, but is in excellent agreement with recent near-IR NICMOS 
observations of the same cluster. 



Subject headings: stars: luminosity function, mass function - globular 
clusters: general - globular clusters: individual: NGC 5139, NGC 6397 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555 
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1. Introduction 

With a total mass of ~ 3.9 10 6 M Q (Meylan 
1987), u;Cen (NGC5139) is the most mas- 
sive globular cluster (GC) in the Galaxy. 
Such a large mass, coupled with the strong 
anisotropy measured in the velocity disper- 
sion (Meylan 1987) implies that, except prob- 
ably in the innermost regions, the cluster has 
not yet reached relaxation through energy 
equipartition and, therefore, the mass func- 
tion (MF) of main sequence (MS) stars should 
still reflect the unperturbed IMF, as no dy- 
namical modifications should have occurred. 
It is, thus, somewhat surprising that Elson 
et al. (1995) obtain for uj Cen a luminosity 
function (LF) which deviates substantially at 
the low mass end from that of other clusters 
with similar metal content. In fact, while 
low metallicity clusters such as NGC6397 
(Paresce, De Marchi, & Romaniello 1995; 
Cool, Piotto, & King 1996), NGC7078 (De 
Marchi & Paresce 1995), NGC7099 (Piotto, 
Cool, & King 1997), NGC 6254 and NGC 6809 
(De Marchi & Paresce 1996), NGC 6656 (De 
Marchi & Paresce 1997) all show a LF that 
reaches a peak at Mj ~ 9 and then drops 
all the way to the detection limit, the results 
published by Elson et al. for uj Cen feature 
a flat LF below Mj ~ 8, with no signs of a 
drop. 

The actual shape of the LF may vary with 
the location inside the cluster because of dy- 
namical evolution, both internal due to two- 
body relaxation and external due to the in- 
teraction with the Galactic tidal field. Yet, 
when the LF is measured near the half-mass 
radius (i.e. within 1-3 times r^ m ) as for all 
the clusters above including to Cen, and bar- 
ring extreme cases of tidal stripping such as 
those resulting from clusters venturing very 
close to the Galactic bulge, one should expect 



only minimal deviations from the original dis- 
tribution of masses as De Marchi & Paresce 
(1996, 1997) have shown observationally, and 
as Richer et al. (1991) had predicted. 

If the results of Elson et al. were con- 
firmed, however, they could modify this pic- 
ture. In fact, because u Cen is both too mas- 
sive to have undergone serious internal dy- 
namical modifications in a Hubble time, and 
is sufficiently far from the Galactic plane and 
bulge not to have experienced major tidal 
stripping (see Dauphole et al. 1996), its LF 
should reflect the IMF much more closely 
than any of the LFs of the other clusters 
listed above. Since clusters with a similar 
metallicity are expected to form with a sim- 
ilar IMF, the discrepancy between the LF of 
uj Cen and that of the other clusters would 
seem to suggest that GCs are born with a 
steeply rising IMF that is substantially modi- 
fied in less massive clusters by dynamical evo- 
lution throughout the clusters lifetime. Thus, 
in these cases the observed MF no longer cor- 
responds to the IMF. 

Such a conclusion would call for a reinter- 
pretation of most of the recent HST-based 
studies of stellar populations in GCs and in 
the field which today all point to a deficiency 
of low-mass stars in the IMF. The implica- 
tions that such a result could bear on the ori- 
gin and nature of dark matter in the Galaxy 
and the star formation mechanism are so im- 
portant that we must be able to exclude that 
the apparent anomaly of the LF of uj Cen re- 
ported by Elson et al. is not due to the 
data reduction process or to the limitation of 
the data themselves. Most of Elson et al.'s 
conclusions are based on the analysis of a 
deep / band image, so that their LF is con- 
structed using little or no color information 
for MS stars. This is somewhat worrisome, 
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as there are many examples of GC LFs based 
on single-filter photometry (e.g. Fahlman et 
al. 1989; Richer et al. 1991) which have later 
proved to be incorrect at the faint end, al- 
though in these cases the error might have 
originated mostly because of crowding. 

To better understand this issue, I have 
used a larger set of WFPC2 images of u> Cen 
now available taken from the HST archive, 
and have subjected these data to exactly the 
same processing that I adopted for all other 
clusters studied so far. In this paper, I de- 
scribe the LF obtained in this way, and com- 
pare it with those of other GCs observed with 
the same instrumentation, as well as with the 
near-IR LF recently obtained with NICMOS 
for the same object (Pulone et al. 1998). 

2. The Data 

The data discussed in this paper have been 
obtained with the WFPC2 on board the HST, 
using the F606W and F814W filters (Biretta 
1996). Two fields have been observed: the 
first, located ~ 7' SW of the cluster center, 
was imaged on 1994 August 29 with a total 
exposure time of 1800 s in the F606W band 
and 2800 s in the F814W filter, while the sec- 
ond, located ~ 4.' 6 SE of the center, was ob- 
served on 1995 June 4 with a total exposure 
time of 3000 s in F606W and 4600 s in F814W, 
respectively. The images of the first field are 
the same as those described in the Elson et 
al. paper (their Field 2), although they erro- 
neously quote for it a distance of 17' 4 from 
the cluster center. The analysis described 
here has been limited only to the PC chips in 
both fields, as the high stellar density makes 
it difficult to obtain accurate photometry for 
faint stars in the WF chips because of their 
larger plate-scale (0"lpix _1 vs 0^043 pix^ 1 ), 
thus increasing the uncertainty on the derived 




Fig. 1. — Top: negative image of a region 
of ~ 30" x 30" located - 7' SW of the center 
of W Cen as seen through the F814W filter of 
the WFPC2 (PC chip) with an exposure of 
1800 s duration (Field 1). Bottom: negative 
image of a region of ~ 30" x30" located ~ 4f 6 
SE of the center of u Cen observed with the 
F814W filter of the WFPC2 (PC chip) with 
an exposure of 4600 s duration (Field 2); the 
box on the lower right corner marks the region 
in which photometry was carried out. 
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LF. The two PC fields used in this paper are 
shown in Figure 1 . 

After standard HST WFPC2 pipeline pro- 
cessing, including frame registration and cos- 
mic ray rejection, stellar fluxes have been 
measured using the core aperture photometry 
technique (De Marchi et al. 1993), with the 
same selection of parameters that I used in 
previous investigations of other clusters with 
the same instrumentation (see Paresce et al. 
1995 and De Marchi & Paresce 1995 for a 
detailed description of our photometric anal- 
ysis and for a discussion of the current un- 
certainties in the reduction and calibration of 
WFPC2 photometric data). 

With a detection threshold conservatively 
set at 5cr above the local average background 
in the PSF peak, I have selected a total of 
1070 stars in Field 1 with well defined fluxes 
in both bands. In Field 2, because of the pres- 
ence of three very bright stars, photometry 
has been restricted to the lower right por- 
tion of the image (see Figure 1), where I have 
found 830 objects satisfying the detection cri- 
terion above in both filters. 

Instrumental magnitudes have been con- 
verted into the WFPC2 ground system by 
using the photometric calibration prescrip- 
tions for the WFPC2 (Holtzman et al. 1995). 
These have been used to construct color mag- 
nitude diagrams (CMDs) extending down to 
m 814 ~ 25, as shown in Figure 2. The inter- 
nal accuracy of the photometry varies in both 
fields from a few hundredths of a magnitude 
for bright, unsaturated stars at m 8 i 4 — 19 
to ~ 0.1 mag at m 8 i 4 ~ 24, but the pho- 
tometry in Field 2 is consistently noisier than 
in Field 1 because of crowding, which makes 
background determination more uncertain in 
a denser environment (the number of stars per 
unit area is more than three times larger in 
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Fig. 2. — Color-magnitude diagrams of the 
stars in Field 1 (left) and Field 2 (right). The 
main sequence in Field 2 is truncated above 
m 8 i4 — 18.5 because of saturation. 

Field 2). This effect can be seen by compar- 
ing the "width" of the MS in both fields and 
noticing that the MS in Field 1 is narrower 
than the other at any given magnitude. The 
MS of uj Cen is truncated at the bright end 
because of saturation, setting in at m 814 ~ 18 
in Field 1 and m 8 i 4 ~ 19 in Field 2 owing to 
the different exposure duration. 

From the CMDs of Figure 2, and adopt- 
ing the 2.5 a clipping criterion described in 
De Marchi & Paresce (1995), I have measured 
the LF of MS stars by counting the number 
of objects in each 0.5 mag bin along the ra 8 i 4 
axis and within ±2.5 times the m 60 6 — m 8 i 4 
color standard deviation around the MS ridge 
line. The contamination due to field stars was 
removed by counting the number of objects 
falling outside of the quoted ±2. 5a belt but 
still within the 0.5 < m 606 — m 814 < 3 color 
range. This number has been scaled by the 
effective color range (2.5 — 5a mag), multi- 
plied by the "width" of the MS belt (i.e. 5a 
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mag) to determine the expected number of 
field stars falling within the cluster MS, and 
finally subtracted from the MS star counts. 
The selection of the above color limits stems 
from the working assumption that field ob- 
jects fill almost uniformly the quoted range. 
The density of objects on both sides of the 
MS in Figure 2 is rather low, in agreement 
with theoretical expectations (Mendez & Van 
Altena 1996). It is, therefore, unlikely that 
we are underestimating the number of stars 
in the MS due to field contamination (see also 
Elson et al. on this issue). 

A reliable determination of the LF from 
the CMDs of Figure 2 requires the assessment 
of the degree of photometric incompleteness 
as a function of the magnitude. I have, there- 
fore, run a number of trials by adding artifi- 
cial stars to the individual frames, and then 
combined and reduced them again with the 
same parameters used in the scientific inves- 
tigation (see De Marchi & Paresce 1995 for 
more details). In this way, I have determined 
that photometry in Field 1 is ~ 100% com- 
plete down to m 814 ~ 22.5, and that, because 
of crowding, the completeness decreases to ~ 
85 % at m 8 i4 ~ 24. Field 2 is more crowded, 
and photometric completeness is ~ 100 % for 
stars brighter than m 814 ~ 20 dropping to 
~ 50 % at m 8 i4 ~ 24. 

3. Discussion and Conclusions 

The LFs measured in this way and cor- 
rected for photometric incompleteness are shown 
in Figure 3, where Field 1 and Field 2 data are 
indicated with boxes and circles, respectively. 
I have applied a vertical offset of 0.1 dex to 
the Field 2 data to account for the different 
stellar density in the two images. In addi- 
tion, I have used different binning for the 
data in Field 1 and Field 2: the size of the 
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Fig. 3. — Luminosity functions of uj Cen mea- 
sured in Field 1 (squares) and Field 2 (crosses) 
using the CMDs of Figure 2. The centers of 
the bins used to count stars along the MS in 
the two fields are displaced by 0.2 mag from 
each other to provide a finer sampling of the 
magnitude range. 

m 8 i4 magnitude bin is the same (0.5 mag) for 
both fields, but the centers of the bins are 
displaced by 0.2 mag from each other. In 
this way, when the LFs are combined in the 
same graph (Figure 3) they provide a finer 
sampling of the magnitude range. Absolute 
magnitudes in Figure 3 have been obtained 
by adopting for co Cen a distance modulus 
of (m — M) = 13.45 and a color excess of 
E(B -V) = 0.15 (Djorgovski 1993), thus im- 
plying (m — M)i = 13.7. The LFs are in ex- 
cellent agreement with each other, within the 
statistical errors, over a range of more than 5 
magnitudes, thus indicating that both fields 
are sampling the same stellar population. 

The shape of our LFs, however, differs sub- 
stantially from that derived by Elson et al., 
in that those presented here reach a peak at 
Mi ~ 9 and then clearly drop all the way 
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to the detection limit, well before photomet- 
ric incompleteness becomes significant and re- 
gardless of whether the field at 4.' 6 or the one 
at 7' from the center is considered. This be- 
havior is perfectly consistent with that ob- 
served in all the other low-metallicity GCs 
studied so far with the WFPC2 (see refer- 
ences in Introduction). To confirm this, in 
Figure 4 I compare the LFs derived here with 
the deepest one today available for a low- 
metallicity GC, namely that of NGC 6397, as 
measured by Paresce et al. (1995) and arbi- 
trarily normalized along the vertical axis. In 
spite of a small difference in the metal content 
of the two clusters ([Fe/H]= —1.59 for cuCen 
and [Fe/H]= -1.91 for NGC 6397; Djorgov- 
ski 1993), their LFs agree very well with each 
other within the observational uncertainties. 
This is not surprising, as both clusters were 
imaged near their half- mass radius (r\ m = 4.' 8 
for c^Cen; Djorgovski 1993), where the LF 
should directly reflect the IMF, as the local 
stellar population is expected to be rather in- 
sensitive to dynamical modifications (Richer 
et al. 1991; De Marchi & Paresce 1997). 

Pulone et al. (1998) have recently observed 
with NICMOS on board the HST a region lo- 
cated ~ 7' SW of the cluster center and have 
derived a LF for MS stars extending down to 
H ~ 22.5 (mieo ^ 26 in the STMAG mag- 
nitude system), corresponding to a mass of 
~ 0.2 M Q . I have converted this LF from the 
H to the / band using the theoretical mass- 
luminosity relation calculated by Baraffe et 
al. (1997) for the metallicity of u Cen, as their 
models of stellar atmospheres have shown to 
reproduce quite accurately the lower portion 
of the MS of low-metallicity GCs (Chabrier 
& Mera 1997). The shape of the /-band LF 
obtained in this way (dashed line in Figure 4) 
is in excellent agreement with those measured 
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Fig. 4. — The luminosity functions of Fig- 
ure 3 are here compared to the LF of the same 
cluster measured with NICMOS (dashed line) 
and converted to the F555W, F814W plane, 
as well as with the I-band LF of NGC 6397. 

in the WFPC2 images over the whole magni- 
tude range spanned by both instruments. 

It remains difficult to understand the dis- 
crepancy between the results presented here 
and those published by Elson et al. for their 
innermost field, unless the flattening that 
they measure at low masses is an artifact due 
to data processing, and particularly to the 
lack of correction for field star and galaxy 
contamination. On the contrary, the color in- 
formation contained in the CMD that I use 
here should make the discrimination between 
cluster and field objects more robust at the 
faint end, where field stars and faint unre- 
solved galaxies begin to outnumber cluster 
objects. Although it would be possible to use 
a model of the Galaxy to estimate and cor- 
rect for field star contamination, the presence 
of unresolved galaxies and QSOs at I ~ 23.5 
and beyond (see e.g. Flynn, Gould & Bahcall 
1996) can cause large uncertainties. 
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In conclusion, my results strengthen the 
scenario suggested by De Marchi & Paresce 
(1996, 1997) that in order for all the low 
metallicity GCs observed so far near the half- 
mass radius to show very similar LFs at low 
masses, they have to be born with the same 
IMF, which is now directly reflected in the LF 
itself with negligible modifications induced by 
dynamical evolution. This proves that it is 
possible to determine the IMF of GC stars, 
provided that a reliable M-L relation is avail- 
able and that the LF is measured near the 
half-mass radius. 

I would like to thank Francesco Paresce 
for many useful discussions on the issues ad- 
dressed in this paper and France Allard and 
Isabelle Baraffe for providing their mass-lumi- 
nosity relations for low-mass stars. 
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